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Experimental  data  on  the  oxidation  kinetics  of  SiC-containing 
diborides  of  Zr  and  Hf  in  the  temperature  regime  of  1473- 
2273  K  are  interpreted  using  a  mechanistic  model.  The  model 
encompasses  counter-current  gas  diffusion  in  the  internal  SiC 
depleted  zone,  oxygen  permeation  through  borosilicate  glass 
channels  in  the  oxide  scale,  and  boundary  layer  evaporation  at 
the  surface.  The  model  uses  available  viscosity,  thermodynamic 
and  kinetic  data  for  boria,  silica,  and  borosilicate  glasses,  and 
a  logarithmic  mean  approximation  for  compositional  varia¬ 
tions.  The  internal  depletion  region  of  SiC  is  modeled  with 
CO/CO2  counter  diffusion  as  the  oxygen  transport  mechanism. 
Data  reported  for  pure  SiC  in  air/oxygen,  for  ZrB2  containing 
varying  volume  fractions  of  SiC,  and  for  SiC-HfB2  ultra-high 
temperature  ceramics  (UHTCs)  by  different  investigations  were 
compared  with  quantitative  predictions  of  the  model.  The 
model  is  found  to  provide  good  correspondence  with  labora- 
tory-furnace-based  experimental  data  for  weight  gain,  scale 
thicknesses,  and  depletion  layer  thicknesses.  Experimental  data 
obtained  from  arc-jet  tests  at  high  enthalpies  are  found  to  fall 
well  outside  the  model  predictions,  whereas  lower  enthalpy  data 
were  closer  to  model  predictions,  suggesting  a  transition  in 
mechanism  in  the  arc-jet  environment. 


I.  Introduction 

Refractory  diborides,  especially  ZrB2  and  HfB2,  with 
SiC  additions  are  being  studied  with  great  interest 
because  of  their  high  thermal  conductivity,  high  melting 
point,  and  moderate  resistance  to  environmental  degradation 
in  air  at  very  high  temperatures.^”^  This  combination  of 
properties  of  these  two-phase  composites  makes  them  prom¬ 
ising  for  use  in  hypersonic  vehicles  as  leading  edge  compo¬ 
nents,  which  are  subject  to  high  heat  fluxes.^  Several 
additions  and  compositional  refinements  are  being  tried  to 
enhance  their  oxidation  resistance,  but  SiC-containing  ZrB2 
or  HfB2  are  the  most  studied,  with  sufficient  data  to  enable  a 
modeling-based  analysis  and  understanding  of  their  oxidation 
behavior.^'** 

A  brief  summary  of  experimental  observations  is  as  fol¬ 
lows. The  oxidation  behavior  is  dominated  by  parabolic 
kinetics  despite  formation  of  a  porous  oxide  of  the  refractory 
metal.  This  protective  behavior  is  attributed  to  the  borosili- 
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cate  glass  that  fills  the  pores  of  the  porous  oxide  scale  and 
often  also  results  in  a  glassy  external  layer,  as  seen  from 
Fig  1(a)  (courtesy  Prof.  E  Opila,  University  of  Virgina).  In  a 
furnace  under  static  air,  a  continuous  external  glassy  layer  is 
formed,  but  at  high  temperatures,  it  begins  to  flow  due  to  a 
decrease  in  viscosity  giving  rise  to  some  scatter  in  experimen¬ 
tal  data.  The  underlying  porous  oxide  is  almost  always  filled 
with  the  borosilicate  glass.  Underneath  this  oxide  layer, 
some  investigations  have  reported  a  depleted  region  where 
SiC  is  absent  but  the  diboride  is  still  intact  as,  for  example, 
seen  in  Fig.  1(a).  The  depleted  region  is  absent  at  low  tem¬ 
peratures,  but  has  been  reported  in  samples  exposed  to  high- 
temperature  air.  However,  there  is  an  inconsistency  in  the 
reported  data  with  some  studies  not  finding  any  depleted 
layers  under  the  same  conditions  where  others  have  observed 
them.^’^^’^"^’^^’^^  In  some  work,  a  partially  depleted  region 
consisting  of  SiC  and  a  Si-O-C  phase,  is  reported.  The  high¬ 
est  temperatures  are  achieved  under  arc-jet  conditions,  where 
typically  very  thick  oxide  layers  with  a  very  thin  or  no  exter¬ 
nal  glassy  layer  are  reported  for  very  short  times.  A  depletion 
zone  has  also  been  observed  under  these  conditions.  Some 
works  have  suggested  that  convection  currents  can  be  signifi¬ 
cant  during  oxidation.  To  obtain  a  quantitative  model 
that  interprets  these  data  was  the  objective  of  this  work. 

In  an  initial  study,  we  reported  on  a  model  for  the  oxida¬ 
tion  of  monolithic  diborides  of  Zr  and  Hf  in  the  temperature 
region  of  1073-2673  In  the  present  work,  we  add  the 

complexities  arising  from  the  presence  of  SiC  and  its  oxida¬ 
tion  products,  silica  and  CO.  We  retain  all  the  essential  ele¬ 
ments  of  the  prior  model,  including  volumetric  effects  from 
phase  change  of  the  oxide,  and  transitions  in  mechanistic 
regimes  as  temperature  is  increased. As  in  the  prior  model, 
the  refractory  oxide  in  the  scale  is  assumed  to  be  imperme¬ 
able  to  oxygen  due  to  the  known  low  electronic  conduc¬ 
tivity.^^  The  effects  of  the  presence  of  SiC  are  the  formation 
of  a  borosilicate  glass  instead  of  a  boria  glass  with  all  its 
consequences,  the  formation  and  transport  of  gaseous  CO 
and  SiO,  and  the  evaporation  of  SiO  and  Si02  at  the  surface. 
The  properties  of  silica  and  boria  are  known  but  modeling 
their  variation  with  composition  needs  interpolation.  The 
vapor  pressures  of  B2O3  and  Si02  are  known  as  a  function 
of  temperature,  and  that  of  SiO  is  known  as  a  function  of 
temperature  and  oxygen  partial  pressure. 

The  model  is  able  to  interpret  most  of  the  experimental 
data  reported  on  ultra-high  temperature  ceramics  (UHTCs). 
Weight  gain,  oxide  thickness,  external  glass  thickness,  and 
internal  depletion  layer  thickness  for  both  SiC-ZrB2  and  SiC 
-HfB2  are  found  to  be  in  reasonable  agreement  with  experi¬ 
mental  data,  but  with  some  exceptions  mostly  related  to 
arc-jet  test  data.  The  equations  that  constitute  the  model  are 
presented  first;  this  is  followed  by  a  section  where  the  predic¬ 
tions  are  compared  with  the  available  experimental  data. 
The  final  section  discusses  the  merits  and  limitations  of  the 
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Fig.  1.  (a)  SEM  image  of  the  microstructure  of  oxidation  scale  formed  on  a  ZrB2-SiC  sample  (Courtesy:  E  Opila,  Univ.  Virginia)  (b,  c) 

Schematic  sketches  of  the  oxidation  products  and  morphology  assumed  in  the  model.  At  lower  temperatures  (b)  external  glassy  scale  forms, 
whereas  at  higher  temperatures  or  high  ambient  flow  (c),  the  glassy  scale  recedes  inward  due  to  evaporative  loss  of  Si02  and  B2O3. 


model,  including  possible  reasons  for  discrepancies  and 
suggestions  for  future  work. 

II.  The  Model 

(1)  Model  Framework  and  Key  Assumptions 
The  morphology  of  the  oxidation  scale  assumed  in  the  model 
was  derived  from  published  cross-sectional  micro  structures  of 
oxidized  The  key  elements  of  the  model 

are  shown  in  Fig.  1.  Figure  1(a)  shows  a  sample  micro¬ 
structure  of  the  scale  formed  on  a  ZrB2-SiC  material.  The 
illustration,  in  Figs.  1(b)  and  (c),  details  the  assumed  mor¬ 
phology  of  the  oxidation  products,  representing  steady-state 
conditions.  The  substrate  is  a  composite  of  SiC  and  MeB2 
(Me  =  Zr,  Hf)  with  being  the  volume  fraction  of  SiC.  The 
oxidation  product,  viz.  the  scale,  consists  of  disconnected  (in 
cross-section)  Me02  grains  with  a  continuous  porous  (void 
of  solid)  region  that  is  filled  with  liquid  borosilicate. 
Throughout  the  manuscript,  pore  refers  to  a  region  void  of 
Me02  but  filled  with  glass,  except  at  very  high  temperatures 
when  glass  evaporates  leaving  pores  behind  [Fig.  1(c)].  The 
Me02  grains  are  taken  to  be  impermeable  to  oxygen,  as 
explained  in  our  prior  work.^^  Molecular  oxygen  from  the 
ambient  dissolves  in  borosilicate  as  molecular  O2,  diffuses 
across  the  external  glassy  phase,  and  permeates  through  the 


porous  channels  of  glass,  eventually  reaching  the  interface  il, 
where  it  reacts  with  the  diboride.  A  portion  of  the  oxygen 
flux,  assumed  to  be  proportional  to  areal  (=  volume)  fraction 
of  SiC,  is  transported  within  the  depleted  region  (void  of 
solid  but  filled  with  gases)  to  oxidize  SiC  through  a  medium 
of  a  gas  mixture  of  CO  and  CO2,  similar  to  the  model  by 
Holcomb  and  St.  Pierre  for  HfC  oxidation.^^  The  gaseous 
reaction  product,  CO,  of  the  SiC  oxidation  is  assumed  to 
either  diffuse  or  bubble  through  the  glassy  scale  (as  observed 
in  ref.  17-19)  and  is  assumed  not  to  be  rate  limiting,  consis¬ 
tent  with  prior  assumptions  for  oxidation  of  SiC.^^  Fig¬ 
ure  1(b)  depicts  the  expected  scenario  at  higher  temperatures 
(likely  above  2000  K)  where  evaporation  of  B2O3  and  Si02 
from  the  external  surface  is  sufficiently  high  that  an  external 
glassy  layer  cannot  be  supported  and  the  glassy  layer  recedes 
inwards.  In  all  regimes,  the  diffusivity  of  gases  in  a  multi- 
component  gas  mixture  and  the  effect  of  Knudsen  diffusion 
[when  an  unfilled  pore  was  present  as  in  Fig.  1(b)]  were  mod¬ 
eled  as  detailed  in  prior  work.^^ 

The  volume  fraction  of  SiC  in  the  substrate  is  /s,  and  the 
volume  fraction  of  porosity  (filled  with  glass)  within  the  Me02 
region  is  /p.  As  in  prior  work,  the  pore  fraction, /p,  (filled  with 
glass)  is  taken  to  change  as  the  temperature  crosses  the  mono¬ 
clinic-tetragonal  phase  transformation  temperature,  Ttrans- 
The  parameters  /§  and  /p  are  related  to  the  volume  fraction. 
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/Me02  5  of  Me02  in  the  glass  +  Me02  region  (2-3),  and  fraction 
of  glass, /g,  (borosilicate)  in  the  scale  as  follows. 


/p  =/ifor  r<rtrans  and /2  for  T  >  Ttrans 

/g  =/s  +/p(l  -/s)  (1) 

/MeO.  =  (1  -/p)(l  -/s)  =  1  -/g 

As  shown  in  Fig.  1,  (1  —  /s)  is  the  volume  fraction  of 
MeB2  and  (1  —  ff)  is  the  volume  fraction  of  Me02  within 
the  Me02  +  glass  layer  (region  2-3).  Thus,  when  no  SiC  is 
present  (monolithic  MeB2),  /g  =  /p.  In  the  presence  of  SiC, 
the  glass  volume  fraction  is  increased  by  f. 

The  B2O3  activity,  of  the  borosilicate  scale  at  the 

MeB2-Me02  interface,  denoted  by  the  subscript  /2,^is  derived 
in  terms  of  the  recession  rate  of  the  MeB2  phase,  — and 
that  of  the  SiC  phase,  Accounting  for  their  phase  frac¬ 
tions  and  their  molar  volumes,  FmcBs  ,  Vsic  of  MeB2  and 
SiC,  and  assuming  an  ideal  mixture,  we  obtain: 
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At  steady  state,  the  fluxes  of  the  gaseous  species  within  the 
depleted  zone  (between  /I  and  /2)  must  be  decided  and  be 
consistent  with  the  equilibrium  constants  of  Eqs.  (4),  and  the 
oxidation  rate  is  limited  by  incoming  oxygen  flux  at  /2, 
J02-32  times  the  area  fraction,  of  SiC  at  the  interface  i2. 
Thus,  for  the  depleted  zone,  we  obtain: 


l-^co-nl  =  2  l•^co2-2l|;  |«/sio-i2|  =  2  l•^co2-2l|; 
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Here,  the  activity  is  taken  to  be  directly  proportional  to 
molar  ratio,  assuming  an  ideal  behavior.  This  assumption  is 
fairly  consistent  with  the  work  of  Boike  et  al,  who  reported 
only  a  slight  positive  deviation  from  ideal  behavior  and  the 
experimental  work  of  Kawamoto  et  al,  who  found  no  phase 
separation  in  this  B203-Si02  system.^^’^^ 

The  equilibrium  constant  for  the  oxidation  of  diboride  is 
related  to  the  boria  activity  and  oxygen  partial  pressure  at 
interface  il  and  is  given  by: 

At  /2  ;  MeB2  +  2  ~  lV[e02  +  B203,’ 

^  ^Me02  ^B203-z2  (2) 

^Me  —  ;  ~TJ^  j 

<^MeB2  (^02-/2)  ^ 

where  K  refers  to  equilibrium  constant,  a  the  activity,  and  P 
the  partial  pressure,  whereas  subscripts  refer  to  species  and 
the  interface  i2  shown  in  Fig.  1 . 


(2)  Internal  Depletion  Region  (1-2) 

Due  to  the  low  oxygen  partial  pressure  at  interface  /2,  the 
only  way  SiC  can  oxidize  at  interface  i\  behind  /2,  is  if 
oxygen  was  transported  across  the  void  by  CO/CO2  counter 
diffusion.  Likewise,  a  gaseous  SiO(g)  from  oxidation  of 
SiC,  transports  Si  to  sustain  the  glassy  oxidation  product  at 
/2.  The  oxygen  dissolved  in  the  glass  at  the  interface  i2 
reacts  with  CO(g)  to  form  C02(g)  which  diffuses  across 
region  (1-2)  to  oxidize  SiC  at  interface  /I.  The  SiO(g) 
formed  at  interface  i\  diffuses  back  to  interface  i2  where  it 
is  oxidized  by  the  dissolved  oxygen  in  the  glass.  This  model 
yields  the  following  reactions  with  associated  equilibrium 
constants. 

At/2  ;SiO(g)  +  ^02  (dissolved)  =  Si02; 

2CO(g)  +  O2  (dissolved)  =  2C02(g) 

3 

Combining  :  SiO(g)  +  2CO(g)  +  -O2  (dissolved) 

=  Si02  +  2C02(g) 

At  i\  ;SiC  +  2C02(g)  =  SiO(g)  +  3CO(g) 


In  these  expressions,  J  is  the  flux  in  moles  per  unit  area 
per  unit  time,  D  the  diffusivity,  R  the  universal  gas 
constant,  T  the  absolute  temperature,  and  /12  is  the  length 
of  the  depleted  zone.  The  factor  10^  refers  to  conversion 
from  atm  to  Pascal.  The  first  subscript  used  in  all  the 
quantities  refers  to  the  species.  The  second  subscript  used 
with  partial  pressures  refers  to  the  interface,  while  that 
used  with  diffusivity  and  flux  refers  to  the  regions  between 
interfaces. 

There  are  seven  unknowns  that  describe  region  1-2,  of 
which  the  oxygen  partial  pressure  at  i2  is  given  by  Eq.  (3). 
The  remaining  six  unknowns,  the  partial  pressures  of  the 
three  species  (CO,  CO2,  and  SiO)  at  the  two  interfaces  i\  and 
/2,  are  obtained  by  solving  the  Eqs.  (4)  and  (5).  Thus,  we 
obtain  the  following  expressions  for  these  six  unknowns. 
They  are  presented  in  algorithmic  sequence,  wherein  each 
variable  can  be  calculated  from  a  knowledge  of  all  preceding 
variables. 
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One  of  them,  Pco-/n  requires  numerical  solution  of  the 
implicit  equation.  The  only  unknowns  in  the  above  equation 
set  are  the  oxygen  flux,  I/02-32I,  in  the  region  2-3  (Fig.  1), 
and  the  length  of  the  depleted  zone,  /12.  The  flux  I/02-32I  is 
obtained  from  solving  the  equations  that  apply  to  region 
2-3,  discussed  below. 


(3)  Me02  +  B203-Si02  Glass  Region  (2-3) 

The  equations  that  govern  the  oxygen  flux  in  region  2-3, 
between  interfaces  i2  and  /3,  are  given  below. 


I  T  I TT  Bo2  —  i3  Po2  —  i2^ 

WO2-32  |— tl02-B203-Si02 - -J - -Jg 

/23 

P02-a~P02-i2> 


\JO2-a3  |— no2-B203-Si02(3fl) 

Since, |/o2-fl3  H-^02-32 


Ha 


Po2-b= 


I23  no2  -B2O3  -Si02  {3a)  P  O2  -fl+/g^3flno2  -B2O3  -Si02  P  O2 
/23no2-B203-Si02(3fl)+/g^3flno2-B203-Si02 


(7) 


Here,  11  refers  to  the  permeability  (rate  of  transport  in 
moles  per  unit  area  per  unit  time  per  unit  partial  pressure 
gradient).  The  permeability  in  region  3-a  is  different  from 
that  in  region  2-3,  as  region  3-a  is  likely  to  be  a  boria-defi- 
cient  region  resulting  from  faster  evaporation  of  boria  than 
silica. 

The  unknowns  that  remain  are  the  three  lengths,  /12,  I23, 
and  Ha  of  the  three  regions.  These  are  allowed  to  evolve  with 
time  from  an  initially  small  value  within  a  numerical  simula¬ 
tion  model.  The  evolution  equations  for  H2  and  derived  from 
flux  and  molar  volumes  are  given  below.  The  rate  of  change 
of  I23  is  given  by  the  number  of  moles  of  Me02  formed  per 
unit  area  per  unit  time  multiplied  by  the  molar  volume  of 
Me02.  The  rate  of  formation  of  Me02  is  given  by  the  total 
flux  of  oxygen  in  through  region  3-2,  less  than  that  used  for 
oxidizing  SiO(g)  to  Si02  and  CO(g)  to  CO2,  which  is  3/2 
times  the  SiO(g)  flux  in  region  1-2  as  given  in  Eq.  4.  The 
recession  rate  of  MeB2  and  SiC  are  simply  related  to  the  rate 
of  change  of  I23  through  molar  volumes  and  volume  frac¬ 
tions.  Finally,  the  depleted  zone  length  change  is  given  by 
the  difference  in  recession  between  MeB2  and  SiC.  Thus,  we 
obtain  the  following  equations. 


d/23 

dt 


FMe02(|-^02-32|  -  2  1-^810-12 


dJ?MeB2 

dt 

dRsic  _ 
dt 


d/23  FMeB2/Me02 
dt  FMe02(l-/s) 


|«JsiO-12| 


Fsic 

/s 


1x2=  7?SiC  —  7^MeB2 


2  1 
5/Me02 


(8) 


In  the  above  set  of  equations,  t  is  time,  R  is  recession,  V  is 
molar  volume,  and  J  is  flux  (moles  per  unit  area  per  unit 
time). 


(4)  External  320^-3102  Glass  Region  (3-a) 

The  evolution  with  time  for  the  thickness,  /3a,  of  the  external 
glassy  layer  3-a,  is  given  by  the  rate  of  production  of  boria 
and  silica  less  the  rate  of  loss  by  evaporation  at  the  surface 
and  the  amount  occupied  by  glass  in  region  2-3.  From  the 
literature,  B203(g),  SiO(g),  and  Si02(g)  are  known  to  be  the 
dominant  gaseous  species  in  the  B203-Si02  system  in  the 


temperature  and  oxygen  partial  pressures  of  interest  (see  for 
example,  ref.  7). 


^Ha  _  ( d/23  /Me02 


dt 


Vmq02 


|•JB203-vap  1^ 


|•Jsi02-vap|  I'^SiO-vap  |)  Fsi02  At 


1^6203  +  (|*JsiO-12| 

d/23 


^  species -vap  | 


TJspecies  ^species— vap 


RT 


^bdry 


dt^^ 

^bdry 


(9) 


/  /specimen 

^flluid 

V  ^fluid 

.Pfluid. 

1/6 


species  ) 


In  the  above  equation  set,  the  subscript  ‘vap’  refers  to 
evaporation,  /specimen  is  the  length  of  the  specimen,  r],  p,  and 
V  refer  to  viscosity,  density,  and  velocity  of  the  ambient  fluid. 
The  last  equation  gives  the  boundary  layer  (thickness  6)- 
dependent  evaporation  of  species. 

The  viscosity  of  the  glassy  layer  could  limit  the  maximum 
thickness,  /3a-max,  of  the  external  layer  that  can  be  supported 
under  gravitational  forces.  As  in  the  previous  work,^"^’^^  from 
the  theory  of  falling  liquid  films,^^  under  laminar  flow,  this 
limiting  thickness  is  given  by: 


Ha—mdiX  — 


^  (^B203  -Si02  rB203  -Si02  /spec  )  ^B203  -Si02 


^/^B203-Si02  sin(0) 


1/3 


rB203— Si02 


dHa  1 

FB203-Si02 


(10) 


where  subscript  B203-Si02  refers  to  the  borosilicate  present 
at  the  surface  region  (3-a),  M  refers  to  molecular  weight,  g 
the  acceleration  due  to  gravity,  cp  the  orientation  of  the  sam¬ 
ple  with  respect  to  gravity,  and  FB203-Si02  refers  to  the  rate 
at  which  the  borosilicate  is  added  to  the  surface  in  moles  per 
unit  area  per  unit  time. 

Finally,  the  net  weight  gained,  weight  of  oxygen  con¬ 
sumed,  and  weight  evaporated  can  be  given  as: 


fFg  —  /23(/Me02PMe02  +/g/^g)  +  HaPg 

—  T^Sic/sPsiC  ~  7?MeB2(l  “/s)PMeB2 

Wo,  =  ^  ^McB,  Mo,  +  ^  Rs,c  f-  Mo, 

^  rMeB2  ^  rsiC 

fFevap  —  'y  ^  Tf species  f  Apecies^^  H  Mqq 

(Si02,Si0,B203) 

(11) 


The  variation  in  time  of  the  oxide  scale  thicknesses  and 
the  various  weight  gain/loss  can  be  computed  numerically 
using  Eqs.  (l)-(ll)  using  an  evolutionary  algorithm  which  is 
started  with  a  choice  of,  initially,  arbitrarily  small  value  for 
/12,  I23,  and  Ha- 


(5)  High-  Temperature  Regime 

At  high  temperatures  (likely  above  2000  K)  and/or  at  high 
flow  rates  of  ambient  fluid  (>1  m/s),  where  the  evaporation 
rates  of  Si02(g)  and/or  SiO(g)  are  sufficiently  high,  the  exter¬ 
nal  glassy  layer  will  be  lost  and  the  glassy  liquid  region  will 
recede  into  the  region  2-3,  as  shown  in  Fig.  1(b).  Using  3/  to 
identify  the  location  of  the  receding  glassy  layer,  the  follow¬ 
ing  equations  describe  the  oxygen  flux  balance  which  fixes 
the  partial  pressure  of  oxygen,  Pq^-bi,  at  the  interface  i3i  at 
the  location  3/, 
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j  T-r  P02-i3i  P02-i2  ^ 

J02-3i2  =  Il02-B203-Si02 - - /g 

j  Po2—3ia  1^5  ^02-a  ~  Po2—i3i  /•  c'  7  t 

J02-a3i  =  Mnce  Jo^-a3i  =  J02-3i2 

p  .  .  —  ^^^^02-B203-Si02^02-i2  -  ^^^^Do2-3iaPo2-a  ~  -^^02-6203-8102-^02-^2 
^^^02-B203-Si02  “  10^^^O2-3w  -  ^^02-B203-Si02 


(12) 


The  diffusive  fluxes  of  gaseous  species  (Si02,  B2O3,  and 
SiO)  determine  the  rate  of  recession  of  the  interface  Bi.  Tak¬ 
ing  the  ratio  of  the  depth  of  the  glassy  region  to  the  length 
I23  of  region  2-3  as  q,  the  evolution  equation  for  q  is  given 
as: 


Js 


D. 


species -vap 


species— 3  zfl 

RT 


10^ 


(1  -  q)l23 


dq  1 
dt  I23 


1  3  2 

7^{*^Si0-12f^Si02  +  {Jo2-3i2  -  ;^«/siO-12)7  ^6203 


-/g 


d/23 


•^Si02-vap  f^Si02  *^B203-vap  ^^6203  }  ^ 


(13) 


The  right-hand  side  of  the  evolution  equation  is  the  rate 
of  production  of  B2O3  and  Si02  less  the  evaporative  flux;  the 
last  term  accounts  for  the  increase  in  I23  during  the  time 
increment.  The  weight  gain  of  the  sample  is  given  as: 


=  ^/23(/Me02PMe02  +/g/^g) 

—  RfsPsiC  ~  ^MeB2(l  — /s)PMeB2 

Again,  the  variation  in  time  of  the  oxide  scale  thickness 
and  the  weight  gain/loss  can  be  computed  using  an  evolu¬ 
tionary  algorithm. 


III.  Model  Predictions  and  Validation 

A  numerical  code  written  in  Fortran  was  used  for  all  the 
model  predictions  shown  in  this  work.  The  input  variables 
are  the  temperature-time  history,  environment  parameters 
(total  pressure,  fraction  of  oxygen),  specimen  orientation, 
length,  fluid  velocity,  volume  fraction  of  SiC,  size  of  SiC  par¬ 
ticles.  All  the  thermodynamic  and  kinetic  parameters  taken 
from  the  literature  reside  within  the  code;  these  include  oxy¬ 
gen  permeation  and  viscosity  of  silica  and  boria,  equilibrium 
constants  for  all  the  reactions,  and  vapor  pressures  of  the 
species.  An  infinitesimal  layer  of  oxide  and  a  glassy  layer 
were  assumed  present  at  an  infinitesimal  start  time  and  the 
time  evolution  of  all  the  parameters  were  calculated  as  dic¬ 
tated  by  the  equations.  An  arbitrary  time-temperature  profile 
could  be  thus  simulated  and  the  resulting  thicknesses  of 
oxide,  glassy  layer,  and  depleted  layer  could  be  computed. 
From  these,  the  recession,  the  weight  gain,  weight  of  oxygen 
consumed,  and  weight  of  evaporated  species  could  be  com¬ 
puted. 


(1)  Parameters 

A  list  of  variables  used  in  the  model  is  shown  in  Table  I  with 
a  brief  description.  The  model  uses  data  from  the  literature 
for  all  thermodynamic  quantities,  viz.  equilibrium  constants 
and  vapor  pressures  of  species,  which  are  available  in  the 
compendium  by  Barin.^^ 

(A)  Gas  Dijfusivity:  The  diffusivities  of  gases  in  a  mul¬ 
tigas  solution  were  calculated  using  parameters  given  by 
Svehla,^^  and  the  methods  outlined  in  Ref.  22  The  gas  diffu- 


sivity  applied  to  the  depleted  region  as  well  as  the  boundary 
layer  diffusion  region  at  the  surface  for  evaporation.  The 
Knudsen  effect  on  diffusivity  was  included  as  determined  by 
the  pore  radius,  assumed  to  be  0.5  pm  in  all  predictions 
shown  herein;  this  choice  was  based  on  reported  microstruc¬ 
tures  in  the  literature.  The  Knudsen  effect  only  comes  in 
when  the  pores  dry  out  under  conditions  where  the  evapora¬ 
tion  is  sufficiently  fast  [Fig.  1(c)].  The  Knudsen  effect  is  also 
calculated  for  the  depleted  zone  as  determined  by  SiC  size. 

(B)  Boria  Activity:  The  calculation  of  equilibrium  par¬ 
tial  pressure  of  oxygen  at  the  ZrB2-Zr02  interface  (/2)  and 
the  partial  pressures  of  volatile  species  at  the  outer  surface 
require  the  activity  of  boria  in  the  borosilicate  glass.  The 
activity  isotherm  in  the  boria-silica  system  was  measured  at 
1475  K  and  is  found  to  be  nearly  ideal,  although  there  is  a 
slight  deviation  from  ideality.^^  The  B203-Si02  system  is  also 
known  for  the  lack  of  any  phase  separation. Based  on 
these,  a  unit  activity  coefficient  was  assumed  in  this  work. 

(C)  Glass  Viscosity:  The  temperature-dependent  vis¬ 
cosity  of  boria  was  obtained  from  the  works  of  Eppler  and 
Li  et  al.^^’^^  The  data  on  viscosity  of  silica  were  reviewed  by 
Doremus,^^  and  as  per  his  conclusion,  the  data  of  Urbain 
et  al?^  were  used  for  the  temperature  range  1400°C-2500°C 
and  that  of  Hetherington  et  al?^  for  the  temperature  range 
1400°C-1000°C.  For  compositional  dependence  of  viscosity 
on  boria  content,  a  log-mean  interpolation  scheme  was  used 
consistent  with  semiempirical  models  of  glass  viscosities. 

(D)  Oxygen  Permeability:  The  permeabilities  of  oxy¬ 
gen  in  liquid  boria  and  silica  were  obtained  by  fitting  to  data 
from  several  sources.  The  boria  data  were  obtained  from 
Tokuda  et  alf^  Luthra,"^^  and  Schlichting."^"^  The  oxygen  per¬ 
meability  in  silica  was  obtained  from  the  works  of  Lamkin 
et  alf^  and  Courtright.^^  For  the  compositional  dependence, 
a  log-mean  approximation  was  used  for  interpolation  consis¬ 
tent  with  what  might  be  expected  from  Stokes-Einstein  rela¬ 
tion  as  suggested  by  Karlsdottir  and  Halloran.^^ 

The  boria  concentration  at  the  surface  layer  must  be  lower 
than  that  in  the  interior  oxide  +  glass  region  due  to  evapora¬ 
tive  losses  of  boria  being  higher  than  silica.  The  surface  layer 
concentration  depends  on  the  evaporative  rates,  diffusion  of 
B  in  B-Si02,  and  the  degree  of  convective  mixing.  This  calcu¬ 
lation  is  beyond  the  scope  of  this  work.  Further,  the  authors 
have  found  no  experimental  data  on  this.  Hence,  the  surface 
boria  concentration  was  assumed  to  be  a  constant  fraction  of 
the  interior  (interior  B2O3  concentration  is  -0.72  for  20% 
SiC);  this  was  the  only  loose  parameter  in  the  model  and  val¬ 
ues  of  0.8  and  0.9  (B2O3  concentration  of  -0.58  and  -0.65) 
for  HfB2  and  ZrB2,  respectively,  gave  the  best  correspon¬ 
dence.  The  effective  pore  fractions,  /p,  for  zirconia  and  hafnia 
were  taken  to  be  the  same  as  was  in  ref.  24  for  monolithic 
diboride  oxidation  (0.03  and  0.04  for  Hf02  and  Zr02,  respec¬ 
tively).  The  aliovalent  dopant  concentration,  Cdopant,  in 
Me02  was  taken  to  be  less  than  100  ppm,  which  permits  the 
neglect  of  oxygen  permeation  through  the  Me02  phase  (an 
excellent  assumption)  (see  ref.  22-24).  The  ambient  fluid 
velocity  was  obtained  from  the  literature  when  reported;  it 
varied  from  0.0001  to  150  m/s,  with  the  smaller  values  corre¬ 
sponding  to  static  air,  and  the  higher  values  close  to  arc-jet 
conditions  (behind  the  shock  wave).  The  effect  of  fluid  flow 
on  evaporation  rates  was  included  in  the  model  using  a 
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Table  I.  A  List  of  Symbols  Used  in  This  Work,  with  a  Brief  Description  and  Units 


Symbol 

Units 

Description 

/p 

Effective  volume  fraction  of  pores  in  Me02  that  is  permeable  to  gas 

fs 

Volume  fraction  of  SiC 

fmeOl 

Volume  fraction  of  Me02  in  the  Me02-glass  region  (2-3) 

A 

Volume  fraction  of  glass  region  in  the  Me02-glass  region  (2-3) 

^species- /2 

Activity  of  species  at  interface  il 

J^MeB2 

m 

Recession  of  MeB2 

^SiC 

m 

Recession  of  SiC 

V 

^  species 

m' 

Molar  volumes  of  species 

P 

^  species 

atm 

Partial  pressure  of  species 

•Jspecies-12 

mol/m^-s 

Flux  of  species  from  interface  1-2 

-^species- 12 

m^/s 

Diffusivity  of  species  in  region  1-2 

R 

J/mol-K 

Universal  gas  constant 

no2 —6203—8102 

mol/m-s-atm 

Permeability  coefficient  of  oxygen  in  liquid  borosilicate 

hi 

m 

Depth  of  internal  depletion  {=R^  —  R^o) 

hi 

m 

Thickness  of  zirconia  region  in  the  scale 

ha 

m 

Thickness  of  external  glassy  layer  of  B203-Si02(/) 

q 

m 

Thickness  of  zirconia  scale  over  which  B203-Si02(/)  is  present 

T 

K 

Temperature 

t 

s 

Time 

Mt 

kg/mol 

Molecular  weight  of  species  i 

9i 

kg/m^ 

Density  of  species  i 

Pa-S 

Viscosity  of  species  i 

FB2O3  — Si02 

mol/m^-s 

Rate  of  addition  of  boria  to  the  external  boria  scale 

•Jspecies-vap 

mol/m^-s 

Rate  of  evaporation  of  species  i  at  the  external  surface 

^bdry 

m 

Boundary  layer  thickness  for  surface  evaporation 

^specimen 

m 

Length  of  specimen 

Ffluid 

m/s 

Velocity  of  ambient  fluid 

kg/m^ 

Net  change  in  weight  per  unit  area 

Wo, 

kg/m^ 

Weight  of  O2  consumed  per  unit  area 

w 

evap 

kg/m^ 

Weight  of  evaporated  species  per  unit  area 

boundary  layer  calculation  (as  detailed  in  ref.  22-24),  but  the 
possible  thinning  of  the  external  glassy  layer  from  shear 
forces  was  not  included  in  the  model. 


(2)  Comparison  of  Model  Predictions  with  Experimental 
Data 

(A)  SiC:  There  are  a  lot  of  oxidation  data  for  SiC  in 
the  literature,  and  comparing  the  model  to  these  data  is 
important  as  a  first  check.  As  the  oxygen  permeability  in  sil¬ 
ica  was  used  in  the  model,  the  comparison  of  model  to  data 
must  be  taken  as  validation  of  the  assumption  that  the  kinet¬ 
ics  of  oxidation  is  limited  by  oxygen  permeability  in  the 
glassy  phase,  at  least  for  SiC.  Figure  2  shows  a  plot  that 
compares  the  available  data  for  single  and  polycrystalline 
SiC  in  various  forms  reported  by  different  investigations"^^“^^ 
in  pure  oxygen;  there  is  a  large  scatter  in  the  data.  Figure  2 
includes  data  collected  on  high  purity  SiC  by  Ramberg 
et  al.f^  and  Ogbuji  and  Opila^^  as  well  as  data  obtained  by 
a  Costello  and  Tressler"^^  on  SiC  made  by  different  processing 
methods.  The  model  is  seen  to  fit  fairly  well  with  the  poly¬ 
crystal  data,  and  the  fast  oxidizing  (Si  face)  data  for  the  sin¬ 
gle  crystal. 

(B)  SiC-ZrB2:  A  comprehensive  work  on  the  effect  of 
SiC  content  on  the  oxidation  kinetics  of  UHTCs  was  per¬ 
formed  by  Talmy.^^  The  study  was  conducted  using  a  ther- 
mogravimetric  apparatus  in  flowing  air.  The  weight  gain  of 
samples  exposed  for  2  h  in  air  is  reported.  In  a  more  recent 
work,  Wang  et  alP  have  conducted  a  similar  study  using 
SiC  volume  fractions  of  0,  5,  10,  15,  and  20  vol%  at  temper¬ 
atures  up  tol600°C  for  Ih  in  an  air  furnace.  Figure  3  com¬ 
pares  the  data  with  the  model  prediction  for  weight  gain  as  a 
function  of  SiC  volume  percent.  Figure  3(a)  shows  data  from 
Talmy^"^  and  Fig.  3(b)  shows  data  from  Wang  et  al}^  The 
model  predictions  for  the  sample  weight  gain  and  the  total 
weight  of  oxygen  consumed  are  shown  as  solid  and  dashed 


1  /  Temperature,  K 

Fig.  2.  The  oxidation  kinetics,  expressed  as  parabolic  rate  constant, 
in  oxygen  of  single  crystal  and  polycrystal  SiC  in  various  forms 
reported  by  various  investigators  are  shown  compared  with  the 
model.  The  model  prediction  is  shown  as  a  solid  line,  whereas  the 
experimental  data  are  shown  as  dotted  lines. 


lines.  The  sample  weight  gain  is  the  weight  of  the  sample  and 
oxidation  products  less  the  weight  of  external  glassy  phase 
predicted  to  be  lost  by  evaporation  or  viscous  flow.  The 
weight  of  oxygen  consumed  will  be  the  total  weight  change  if 
all  of  the  oxidation  products  were  retained  on  the  sample. 
Figure  3(c)  includes  time  dependence  data  at  1200°C  from 
Wang  et  al.  The  model  agrees  with  the  trends  very  well. 

The  temperature  dependence  of  weight  gain  and  scale 
thickness  data  on  a  fixed  composition  of  20  vol%SiC-ZrB2 


344 


Journal  of  the  American  Ceramic  Society — Parthasarathy  et  al. 


VoL  95,  No.  1 


100  200 

Time,  min 


10  15  20 

Volume  %  SiC 


25 


Fig.  3.  (a)  The  effect  of  volume  percent  SiC  on  the  oxidation  weight  gain  in  2  h  of  UHTC  (SiC-ZrB2)  samples,  as  measured  by  Talmy^"^  shown 

compared  with  the  prediction  of  the  model  for  three  different  temperatures.  The  solid  lines  are  the  predictions  that  include  the  flow  of  external 
glassy  layer  under  gravity  and  evaporation  of  boria  or  silica.  The  dotted  lines  show  the  predicted  weight  of  oxygen  consumed,  which  is  the 
maximum  weight  gain  that  the  samples  could  have  suffered  from  oxidation  in  the  absence  of  liow/evaporation  of  the  glassy  layer.  In  (b),  the 
model  is  shown  compared  with  data  from  a  different  investigation,  by  Wang  et  which  used  a  1-h  hold  at  temperature.  In  (c),  the  weight 
gain  as  a  function  of  time  for  a  ZrB2-15  vol%SiC  sample  at  1200°C  reported  by  Wang  et  al}^  is  shown  compared  with  the  model. 


in  static  laboratory  air  have  been  reported  by  Carney  et  al}^ 
This  study  took  special  care  to  document  the  effect  of  glass 
flow  in  two  ways.  They  measured  the  weight  changes  of  both 
the  sample  and  the  crucible  onto  which  some  of  the  glass 
had  spread.  They  also  measured  the  variation  in  scale  thick¬ 
ness  with  sample  surface  orientation,  and  found  significant 
scatter.  Further,  they  reported  on  samples  prepared  by  SPS 
and  hot  press.  For  comparison  with  the  model,  the  average 
values  of  all  the  samples  for  a  given  condition  were  used. 
The  experimental  data  points  attributed  to  Opila  and  Halbig 
in  Fig.  4(a)  were  calculated  from  the  parabolic  constants  for 
weight  gain  reported  in  their  work.^^  In  Fig.  4(a),  the  total 
weight  gain  (sample  +  crucible)  is  plotted  along  with  the 
model  predictions,  showing  good  correspondence.  The  plot 
also  shows  that  evaporation  is  significant  above  1873  K. 
Figure  4(b)  shows  the  average  values  for  total  scale  (oxide 
+  glass)  thicknesses  and  oxide  thicknesses.  Figure  4(b) 
includes  scale  thicknesses  data  measured  from  images 
reported  by  Zhang  et  al}^  The  data  of  Carney  were  obtained 
in  an  alumina  furnace  whereas  the  data  of  Zhang  et  al.  were 
obtained  using  a  zirconia  furnace.  A  clear  internal  depletion 
layer  was  not  observed/reported  in  the  experiments.  In 


contrast,  Fahrenholtz  ^  reported  a  depletion  layer  of  10  pm 
after  30  min  at  1773  K  in  a  ZrB2-30%  SiC  sample;  the 
model  predicts  only  0.6  pm  for  these  conditions.  For  the 
20%  SiC-ZrB2  case,  the  model  predicts  a  depletion  zone 
ranging  from  1.8  pm  at  1673  K  to  3  pm  at  1873  K  which  is 
closer  to  observations  by  Carney  et  al.  Figure  5  shows  a 
comparison  of  the  model  with  data  for  the  time  dependence 
of  weight  gain  and  scale  thickness  from  two  sources  of 
data.^’^^  Data  for  hot-pressed  samples  and  SPS-processed 
samples  and  the  sources  of  the  data  are  distinguished  with 
different  symbol  shapes.  Once  again,  there  is  a  reasonable 
correspondence  between  the  data  and  the  model. 

(C)  SiC-HfB2:  Oxidation  data  on  SPS-processed 
high-density  20  vol%SiC-HfB2  samples  have  been  collected 
in  static  air  from  1673  to  2273  K  by  Carney.  In  addition, 
the  statistical  variations  from  batch  to  batch  and  processing 
route  (hot-pressed  versus  SPS)  have  been  studied  for  the 
same  composition  at  1773  and  1873  K  by  Sevener.^^ 
Figure  6(a)  compares  the  data  for  weight  gain  with  the 
model  predictions.  The  data  for  weight  gain  include  separate 
measurements  for  the  sample  (open  symbols)  and  crucible 
+  sample  (filled  symbols).  The  model  predictions  include 
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Fig.  4.  (a)  Oxidation  weight  gain  in  ZrB2-20  vol%SiC  samples, 

measured  in  static  air  after  2.5  h  of  exposure,  by  Carney  et  al}^ 
shown  compared  to  the  model  predictions.  The  data  from  Opila  and 
Halbig^^  are  also  shown.  The  solid  lines  are  model  predictions  for 
weight  gain,  weight  of  oxygen  consumed,  and  weight  evaporated,  (b) 
Average  of  reported  values  for  the  oxide  scale  thickness  (open 
symbols)  and  total  (oxide  +  glass)  scale  thickness  (filled  symbol)  for 
the  same  exposures;  values  measured  on  different  sides  of  the 
samples  were  averaged  for  this  plot.  Values  from  the  work  of  Zhang 
et  which  used  a  zirconia  furnace,  are  also  included.  The  solid 
lines  are  predictions  for  total  scale,  oxide  scale,  and  depletion 
thicknesses.  No  depletion  thickness  was  reported,  but  partial 
depletion  was  noticed  under  FIB/TEM  investigation.^^ 


Time,  min 

Fig.  5.  (a)  Model  compared  with  oxidation  weight  gain  measured 

as  a  function  of  time  during  oxidation  of  ZrB2-20  vol%SiC  samples 
in  static  air  at  1873  K,  by  Carney  et  and  at  1900  and  1600  K 
by  Levine  et  al.^  The  solid  lines  are  model  predictions  for  the  two 
temperatures,  (b)  Model  compared  with  reported  oxide  thickness 
(open),  total  (oxide  +  glass)  thickness  (filled),  and  internal  depletion 
thickness  (dash)  as  marked.  In  both  (a)  and  (b),  the  squares 
represent  Carney’s  hot-pressed  samples,  the  diamonds  Carney’s  SPS 
samples,  and  triangles  data  from  Levine  et  al. 


total  weight  of  oxygen  consumed,  sample  weight  gain,  and 
weight  evaporated.  Evaporation  is  predicted  to  dominate  at 
around  2000  K,  and  the  model  is  consistent  with  the  drop  in 
sample  weight  gain  at  this  temperature.  In  general,  the  model 
captures  the  trends  well  except  at  the  highest  temperature.  In 
Fig.  6(b),  the  scale  thicknesses  measured  are  compared  with 
the  model.  The  micro  structures  of  the  oxidation  product 
were  complex  at  the  highest  temperatures,  and  there  was 
ambiguity  in  the  definition  of  the  depleted  zone.  Thus,  the 
sum  of  oxide  scale  and  depletion  zone  (open  symbols)  was 
used  for  this  comparison.  The  total  scale  thicknesses  (filled 
symbols)  are  also  plotted.  The  correspondence  between  the 
model  and  data  is  seen  to  be  reasonable,  except  for  the  data 
at  the  highest  temperature  reported,  viz.  2173  K,  which  devi¬ 
ates  significantly  from  the  model.  However,  these  data  also 
deviate  significantly  from  the  extrapolation  of  the  data  at 
lower  temperatures.  On  further  examination  of  the  sample, 
significant  contamination  of  alkali  elements  (mainly  Ca)  pos¬ 
sibly  from  the  furnace  or  crucible  (Ca-stabilized  zirconia) 


was  detected.  Enhanced  oxidation  of  SiC  due  to  the  presence 
of  impurities  is  well  known  (for  example,  see  Ramberg  et  alf^) 

Figure  7  compares  the  model  predictions  for  weight  gain 
as  a  function  of  time  for  20  vol%SiC-HfB2  in  furnace  air  at 
1773  and  1873  K,  with  data  reported  by  Sevener.^®  The  data 
include  SPS-processed  material  and  hot-pressed  material  with 
different  initial  SiC  particle  sizes  obtained  from  different  pro¬ 
cessing  routes.  The  data  reported  for  the  combined  weight 
gain  of  sample  and  crucible  were  used.  The  correspondence 
is  once  again  found  to  be  good,  although  there  is  consider¬ 
able  scatter  in  the  data.  The  model  prediction  plotted  in 
Fig.  7  includes  the  total  weight  of  oxygen  consumed  showing 
that  evaporative  loss  becomes  significant  at  1873  K. 

(D)  Arc- Jet  Tests  on  SiC-MeB2:  The  expense  of  con¬ 
ducting  arc-jet  tests  has  limited  the  number  of  investigations 
and  the  extent  of  data  available  on  oxidation  kinetics  during 
these  tests.  However,  data  are  available  for  a  few  SiC-MeB2 
compositions.  Montverde  and  Savino  have  tested  a  hemi¬ 
spherical  sample  of  ZrB2-15  vol%SiC  under  arc-jet  condi¬ 
tions  where  the  surface  temperature  reached  -2193  K  for 
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Fig.  6.  Comparision  of  model  with  experimental  data  on  oxidation 
kinetics  of  HfB2-20  vol%SiC  in  the  high-temperature  regime  (up  to 
2173  K)  in  static  air  showing  weight  gain  in  (a)  and  scale  thicknesses 
in  (b).  All  the  data  are  from  the  work  of  Carney.  The  solid  lines 
are  model  predictions.  In  (a),  both  sample  weight  (open)  and 
sample  +  crucible  (filled)  weights  are  shown,  indicating  flow  is 
significant  above  2000  K.  Similarly  in  (b),  total  scale  thickness 
(filled)  and  oxide  scale  thickness  (open)  are  shown.  Due  to 
experimental  ambiguity  in  demarcation  between  depletion  and  oxide 
scale,  oxide  scale  +  depletion  thickness  are  shown  for  data  and 
model. 


325  min.^^  From  their  metallographic  images  along  with 
EDS  maps,  the  thicknesses  of  the  oxide  layer,  glass  layer, 
and  depletion  layer  were  obtained.  These  are  plotted  along 
with  model  predictions  in  Fig.  8(a).  A  bilinear  approximation 
of  the  measured  temperature-time  profile,  used  as  input  for 
the  model,  is  shown  as  inset.  More  recently,  Savino  et  al 
have  conducted  a  similar  study  at  higher  surface  tempera¬ 
tures  using  the  arc  jet.^^  From  the  reported  thermal  history 
and  calculated  scale  thicknesses  from  the  published  micro¬ 
graphs,  the  plot  shown  in  Fig.  8(b)  was  generated  which 
shows  the  model  predictions  compared  to  the  data. 

Data  on  arc  jet  tested  samples  of  SiC-HfB2  have  been 
reported  on  by  Carney^  ^  and  Gasch  et  al}^  In  particular, 
Carney  reported  on  scale  thicknesses  for  both  furnace- 
exposed  and  arc-jet-exposed  samples  of  the  same  batch  of 
sample  under  the  same  temperature  (1773  K)  and  duration. 
Figure  8(c)  shows  the  model  predictions  compared  with  fur¬ 
nace  data  and  the  arc-jet  data.  Figure  8(d)  shows  the  data 
from  Gasch  et  al.,  who  conducted  two  tests:  one  at  surface 


Fig.  7.  Oxidation  data  on  different  batches  of  hot-pressed  (a  and  c) 
and  SPS-processed  samples  (b)  of  20  vol%SiC-HfB2  composition, 
collected  at  1773  and  1873  K  by  Sevener.^^  The  solid  lines  are  model 
predictions  for  the  weight  gain  and  weight  of  oxygen  consumed;  the 
difference  arises  from  evaporation  of  volatile  oxides. 


temperature  of  1803  K,  and  another  at  1963  K.  The  model 
predictions  for  both  these  temperatures  are  shown. 

In  general,  the  model  underpredicts  the  kinetics  for  experi¬ 
mental  data  on  samples  exposed  to  arc  jet.  The  model  is 
much  closer  at  lower  temperatures.  The  data  from  furnace- 
exposed  samples  show  slower  kinetics  and  are  in  reasonable 
correspondence  with  the  model. 

IV.  Discussion 
(1)  Model  Strength  and  Weaknesses 

A  model  has  been  presented  to  interpret  the  oxidation  kinet¬ 
ics  of  SiC-containing  refractory  metal  diborides.  The  model 
includes  the  effect  of  viscous  flow  of  the  outer  glassy  scale 
and  evaporative  losses  of  volatile  species  under  a  boundary 
layer  condition  established  by  ambient  fluid  flow  velocity. 
The  model  predicts  the  total  weight  of  oxygen  consumed, 
sample  weight  gain,  weight  of  evaporated  species,  oxide 
thickness,  glass  thickness,  depletion  layer  thickness,  and 
recession  of  the  substrate.  The  model  uses  as  input  parame¬ 
ters,  the  exposure  temperature  and  time  or  a  thermal  profile, 
ambient  oxygen  partial  pressure,  gas  chemistry,  and  fluid 
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Fig.  8.  Results  from  arc-jet  testing  done  on  ZrB2-SiC  (a  and  b)  and  on  HfB2-SiC  (c  and  d)  are  shown  along  with  model  predictions  for  the 
reported  thermal  history  used  in  the  experiment  (shown  as  inset  in  a  and  b)  (data  from  ref.  11,12,55,56).  In  general,  the  model  is  seen  to 
underpredict  the  oxidation  kinetics,  with  the  discrepancy  increasing  with  temperature  of  exposure.  In  one  case,  shown  in  (c),  a  furnace  oxidation 
and  arc -jet  run  were  conducted  for  the  same  temperature  and  time.  The  model  is  found  to  be  much  closer  to  the  furnace-based  data  and  well 
below  the  arc-jet  data. 


velocity.  The  material  parameters  for  input  are  the  size  and 
volume  fraction  of  SiC,  sample  length  along  the  fluid  flow 
direction.  The  model  calculates  the  activity  of  boria  in  the 
oxide  +  glass  region,  but  the  activity  of  boria  at  the  external 
glassy  layer  is  assumed  to  be  a  fraction  of  the  inner  region. 
This  was  the  only  fit  parameter.  The  model  used  a  CO/CO2 
counter  diffusion  mechanism  for  oxygen  transport  within  the 
depleted  region. 

In  Fig.  2,  the  model  is  seen  to  be  in  agreement  with  the 
kinetics  for  oxidation  of  pure  SiC  in  dry  oxygen,  although 
there  is  significant  scatter  in  the  data.  This  agreement  shows 
that  oxidation  kinetics  of  SiC  is  limited  by  diffusion  of  molec¬ 
ular  oxygen  across  the  silica  scale,  assumed  in  the  model.  It 
further  lends  credibility  to  the  assumption  that  the  escape  of 
CO,  one  of  the  oxidation  products,  is  not  rate  limiting. 

From  Fig.  3,  the  effect  of  volume  fraction  of  SiC  on  the 
kinetics  of  ZrB2-based  UHTCs  is  seen  to  be  in  agreement 
with  the  model  for  the  two  different  sources  of  experimental 
data.  The  model  predictions  for  the  total  oxygen  consumed 
are  shown  by  dotted  lines,  and  the  predicted  sample  weight 
gain  by  solid  lines.  The  loss  of  material  from  the  surface  due 
to  evaporation  and  fluid  flow  accounts  for  the  difference 
between  the  two.  Figure  3(c)  confirms  that  the  model  cap¬ 
tures  the  parabolic  time  dependence  quantitatively. 

In  Fig.  4(a),  the  experimental  data  on  weight  gain  agree 
well  with  the  model,  and  the  model  predicts  accelerated  evap¬ 
oration  at  temperatures  above  1873  K.  In  Fig.  4(b),  the  mea¬ 
sured  scale  thicknesses  are  in  agreement  with  model  within 
the  bounds  of  scatter  in  the  experimental  data.  The  scatter 
arises  from  the  variations  in  the  flow  of  externally  glassy 
layer,  as  detailed  in  Carney.  The  data  of  Zhang  et  al}"^ 
were  obtained  using  a  zirconia  furnace,  where  the  possibility 


of  contamination  resulting  in  accelerated  oxidation  cannot  be 
ruled  out.  The  model  predicts  the  depletion  layer  to  be  2  pm 
at  1673  K  rising  to  only  3  pm  at  1873  K.  Figure  5  shows 
that  the  evolution  in  time  of  weight  gain  and  scale  thicknesses 
agree  well  with  the  model.  The  data  of  Carney at  1873  K 
and  the  data  of  Levine  et  al.^  at  1600  and  1900  K  are  used 
to  compare  with  the  model.  The  two  symbols  for  data  of 
Carney  refer  to  samples  processed  by  SPS  and  hot  pressing. 
The  data  of  Levine  et  al  fall  slightly  below  these.  The  deple¬ 
tion  depths  reported  by  Levine  et  al.  at  1900  K  are  higher 
than  the  model  prediction,  while  depletion  was  not  observed 
by  Carney  at  1873  K. 

The  model  predicts  a  finite  but  small  depleted  zone.  The 
experimental  reports  show  large  variations.  A  10- pm  thick 
depletion  zone  was  reported  bv  Fahrenholtz  at  1773  K  in 
30  min,  whereas  Zhang  et  al^"^  report  no  depleted  zone 
below  1873  K.  Similarly,  no  clear  depleted  zone  was 
observed  and  reported  by  Carney  after  2.5  h  at  1673- 
1873  K.  However,  extensive  EDS  and  FIB/TEM  analyses  did 
reveal  a  partially  oxidized  internal  zone  which  was  2  pm  at 
1673  K  and  14  pm  at  1873  K.  As  the  actual  fraction  of  SiC 
oxidized  in  this  partially  oxidized  zone  was  not  measured,  it 
is  difficult  to  compare  with  this  model.  In  summary,  some 
ambiguity  remains  with  respect  to  internal  depleted  zone  for¬ 
mation  in  these  materials. 

In  Fig.  6,  the  HfB2-SiC  data  of  Carney^  ^  are  rationalized 
well  by  the  model  up  to  2073  K,  but  there  is  a  significant 
increase  in  measured  oxidation  compared  to  the  predictions 
at  2173  K.  The  reasons  for  this  remain  unknown,  but  per¬ 
haps  contamination  from  the  stabilized  zirconia  heating  ele¬ 
ments  may  have  caused  this  discrepancy  because  significant 
Ca  was  detected  in  EDS  of  the  scales.  The  scales  were  found 
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to  have  a  complex  morphology  at  temperatures  of  1973  K 
and  above.  The  expected  depletion  layer  was  not  observed; 
several  different  regions  as  detailed  in  Carney  were  noted. 
Thus,  the  plot  only  shows  the  combined  thickness  of  oxide 
and  depletion  layer,  and  the  total  scale  thickness  for  the 
experimental  data.  Within  experimental  scatter,  in  Fig.  7,  the 
model  captures  the  evolution  in  time  of  weight  gain  at  1773 
and  1873  K,  as  reported  by  Sevener.^^ 

Figure  8  reveals  the  key  inadequacy  of  the  model.  The  oxi¬ 
dation  kinetics  during  arc-jet  testing  were  measured  and 
reported  by  Montverde  and  Savino  for  ZrB2-15SiC  under  two 
different  heat  flux  conditions. The  actual  thermal  profiles 
reported  were  used  in  the  model.  The  model  fails  to  capture 
even  the  correct  order  of  magnitude  of  scale  thicknesses  at 
high  heat  fluxes,  while  it  is  reasonable  for  low  heat  flux  condi¬ 
tions.  The  internal  depletion  is  considerably  underpredicted. 
Similarly  for  HfB2-20SiC,  the  model  underpredicts  the  kinet¬ 
ics  by  a  significant  margin.  In  one  reported  work,  both  furnace 
oxidation  and  arc-jet  oxidation  were  conducted  at  the  same 
temperature  and  time,  viz.  1773  K  for  4  min.  Figure  8(d) 
shows  that  once  again,  the  model  captures  arc-jet  data  better 
at  lower  temperatures,  but  not  at  higher  temperatures  using 
data  from  Gasch  et  al.  on  HfB2-20SiC.^^  The  model  predic¬ 
tion  is  close  for  the  furnace  oxidation  data,  but  far  below  the 
arc-jet  data.  Clearly,  aspects  of  oxidation  kinetics  that  take 
place  during  arc-jet  conditions  are  not  captured  in  the  model. 


(2)  Future  Directions 

The  model  is  clearly  lacking  in  capturing  high-temperature 
oxidation  kinetics  under  arc-jet  conditions.  Several  possible 
mechanisms  can  be  suggested  that  require  further  analysis 
and  experimentation.  The  dissociation  of  gases  into  mon¬ 
atomic  oxygen  and  nitrogen  under  arc-jet  conditions  is  well- 
known.(for  example,  see  ref.  57)  However,  the  oxidation 
mechanisms  under  such  conditions  have  only  been  studied  by 
a  few.^^’^^  These  studies  show  that  oxidation  is  indeed  faster 
under  dissociated  oxygen  conditions;  however,  the  reasons 
for  this  have  not  been  established.  The  recombination  of  the 
monatomic  gases  at  the  material  surface  can  release  energy 
that  is  sufficient  to  raise  the  material  surface  temperature.  In 
the  furnace  oxidation  tests  using  monatomic  oxygen,  the 
material  surface  temperature  could  not  be  measured.  How¬ 
ever,  during  arc-jet  testing,  the  material  temperature  was 
measured,  and  the  model  is  unable  to  capture  the  oxidation 
kinetics  with  these  measured  temperatures,  indicating  that 
some  other  mechanism  might  play  a  role.  The  possibility  of 
enhanced  oxidation  kinetics  from  permeation  of  oxygen  in 
monatomic  form  has  been  suggested,  but  no  direct  experi¬ 
mental  evidence  exists.  A  modeling  work  by  Li  et  al.  based 
on  first  principles  showed  that  the  dissolved  concentration  of 
monatomic  oxygen  in  borosilicate  glasses  can  be  higher  than 
that  of  diatomic  oxygen,  at  temperatures  above  1773 
Given  that  the  arc-jet  test  data  are  closer  to  the  model  at 
lower  temperatures  (below  1900  K)  and  deviate  significantly 
at  higher  temperatures,  the  enhanced  permeation  by  mon¬ 
atomic  oxygen  in  borosilicate  glasses  at  higher  temperatures 
is  a  possibility.  Monatomic  oxygen  would  also  change  the 
thermodynamics  of  oxidation  at  the  ZrB^/oxide  interface. 
Theoretical  and  experimental  investigations  that  could  con¬ 
firm  this  hypothesis  are  suggested  for  future  work.  It  is  also 
suggested  that  the  possibility  of  active  oxidation  of  SiC  from 
the  surface  by  monatomic  oxygen  be  investigated.  The  effect 
of  mechanical  shear  at  high  velocities  has  been  neglected  in 
this  work;  however,  the  model  predicts  evaporation  to  be 
high  enough  to  predict  the  absence  of  external  glass  under 
arc-jet  conditions. 

Another  limitation  of  the  model  is  its  inability  to  predict 
the  depletion  layer  thicknesses  in  all  the  reported  works. 
However,  there  is  considerable  disagreement  between  differ¬ 
ent  works  on  the  extent  and  even  existence  of  the  depletion 
layer  under  identical  conditions  for  similar  or  identical  com¬ 


positions.  The  model  shows  clearly  that  counter  diffusion  via 
the  CO/CO2  mixture  as  oxygen  transport  medium  is  tenable 
and  can  lead  to  a  depleted  region.  However,  the  extent  is  not 
as  high  as  in  some  of  the  reports,  while  it  is  consistent  with 
others.  Some  uncontrolled  variable,  such  as  humidity  or  sam¬ 
ple  porosity,  within  the  experiment  might  be  responsible  for 
this.  Future  studies  that  clarify  this  might  be  useful.  Finally, 
the  effect  of  mechanical  forces  that  induce  shear  of  the  exter¬ 
nal  glass  by  viscous  flow  has  been  ignored  in  this  work  and 
may  be  worth  examining  in  future  models. 

V.  Summary 

A  mechanistic  model  has  been  built  to  interpret  experimental 
data  on  the  oxidation  kinetics  of  SiC-containing  diborides  of 
Zr  and  Hf  in  the  temperature  regime  of  1473-2473  K.  The 
model  uses  available  viscosity,  thermodynamic  and  kinetic 
data  on  borosilicate  glasses,  and  uses  a  logarithmic  mean 
approximation  for  compositional  variations.  The  internal 
depletion  region  of  SiC  is  modeled  using  counter  diffusion  in 
CO/CO2  as  oxygen  transport  mechanism.  Data  obtained 
after  furnace  oxidation  of  pure  SiC  in  oxygen,  SiC-contain¬ 
ing  ZrB2  with  varying  volume  fractions,  and  SiC-HfB2 
UHTCs  were  found  to  agree  reasonably  well  with  the  model. 
The  model  predictions  include  weight  of  oxygen  consumed, 
sample  weight  gain,  external  glassy  layer  thickness,  inner 
oxide  scale  thickness,  recession,  and  depletion  layer  thick¬ 
ness.  Some  depletion  layer  thicknesses  and  data  obtained 
using  arc-jet  tests  fall  well  outside  the  model  predictions, 
indicating  that  some  aspects  of  arc-jet  conditions  are  not 
captured  in  this  model. 
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